Equation of state of looped DNA 



Igor M. Kulic 1 , Herve Mohrbach 2 , Rochish Thaokar 3 , and Helmut Schiessel 2 
1 Department of Physics and Astronomy, University of Pennsylvania Philadelphia, PA 19104, USA 
2 Instituut-Lorentz, Universiteit Leiden, Postbus 9506, 2300 RA Leiden, The Netherlands 
3 Department of Chemical Engineering, Indian Institute of Technology, Bombay, Mumbai, 400076, India 

(Dated: February 9, 2008) 

We calculate the equation of state of DNA under tension for the case that the DNA features loops. 
Such loops occur transiently during DNA condensation in the presence of multivalent ions or sliding 
cationic protein linkers. The force-extension relation of such looped DNA modelled as a wormlike 
chain is calculated via path integration in the semiclassical limit. This allows us to determine 
rigorously the high stretching asymptotics. Notably the functional form of the force-extension curve 
resembles that of straight DNA, yet with a strongly renormalized apparent persistence length. That 
means that the experimentally extracted single molecule elasticity does not necessarily reflect the 
bare DNA stiffness only, but can also contain additional contributions that depend on the overall 
chain conformation and length. 

PACS numbers: 87.15.La, 82.37.Rs, 05.40.-a 



I. INTRODUCTION 



The DNA double helix is the molecule that encodes the genetic information in living cells. In addition to carrying 
the genome, DNA has also specific physical properties that are essential for its biological functions. Its mechanical 
properties are exploited by the protein machinery for the transcription, replication, repair and packaging of DNA 0]. 
During the last decade it has been possible to manipulate single DNA molecules to determine its elastic properties 
under different physical conditions 2]. In these experiments, the extension of single molecule versus an applied 
stretching force is measured by a variety of technics including magnetic beads 0, Q , optical traps 0, , micro- 
needles pj, hydrodynamic flow [8| and AFM |j. The studies also made it possible to better understand mechanical 
interactions between DNA and proteins 01 ■ 

The most appealing theoretical description of the DNA molecule is the wormlike chain (WLC) model which is a 
coarse-grained model of DNA with a single parameter, the persistence length l p , characterizing the chain stiffness. 
Originating back to the first half of the last century [TJ it gained renewed interest after the semiflexible nature of 
DNA and other (bio)polymers became clear 0] and it is now indispensable for the theoretical understanding of many 
single molecule experiments. In many cases stiff polymers show a characteristic force-extension behavior that can be 
well understood in terms of the WLC as being the result of entropic fluctuations of the chain that - with increasing 
tension - become suppressed at shorter and shorter wave lengths. Measuring the force-extension characteristics of 
such a chain allows to extracts its overall contour length as well as its persistence length. 

The DNA in living cells is rarely found in its straight "naked" state; rather an overwhelming fraction of DNA is 
strongly configurationally constrained by binding proteins causing loops, bends and wraps. In particular, protein 
complexes forming loops are essential for biological process like distant gene expression or DNA packaging The 
formation of loops in a DNA molecule under tension has been the subject of the theoretical investigation in Ref. . 
Also, single molecule stretching experiments on DNA condensed with multivalent counterions performed by several 
groups 0, ^| might bear loops or related structures like DNA toroids 0]. While the statistical mechanics of 
unconstrained DNA under tension is well described by the WLC 0, the presence of topological constraints like 
supercoiling 0, 0] and entanglements , or geometrical constraints like protein induced kinks and bends [2(], l2l| 
renders analytical results more difficult. 

In this paper, we expand the repertoire of analytically solvable "equations of state" by deriving the force extension 
relation for a DNA with a sliding loop as depicted in Fig. ^ The computation is performed by evaluating quadratic 
fluctuations around the looped solution - a non-constant saddle-point of the DNA elastic energy. The method is 
essentially analogous to the semiclassical treatment of tunneling amplitudes in quantum mechanics and instantons 
in quantum field theory [2^]. The equation of state of looped DNA that we present here can be considered as a 
paradigmatic model case for stretching DNA with a non-trivial ground state. Having understood the physics of the 
looped DNA it is straightforward to extend the analysis also to other cases where the overall DNA conformation 
is far from being straight. The calculation presented in the present work has been sketched in a previous paper 
|23|| (for the high-force limit) together with some interesting experimental situations, namely rigid protein-induced 
kinks and anchoring deflections in AFM stretching of semiflexible polymers. Expressions relating the force-extension 
measurements to the underlying kink/boundary deflection geometry were also provided in Ref. |23| and applied to 
the case of the GalR-loop complex [24| . 
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FIG. 1: Various examples of stable loops in DNA under tension: (a) Applied torque M at the ends, (b) DNA adsorbed on a 
surface, (c) DNA surrounded by a dense solution of infinitely long DNAs. Unfolding of the loop goes hand in hand with an 
energetically costly transient "cavity" creation (in the pink region), (d) DNA in a strong magnetic field H perpendicular to the 
applied force, (e) DNA looped by a freely sliding linker ligand ("weakly condensed" DNA). 

It is not the purpose of this paper to analyze concrete experimental setups, but to give a more detailed and 
comprehensive description of our computation. This paper is also more general as we employ our method also to 
determine the equation of state of sliding loop in the small force regime where the semi-classical approximation is 
valid for DNA with length smaller than the persistence length whereas for longer DNA our computation is only valid 
for large forces. In that strong force regime, it will be shown that the presence of the loop modifies the the elastic 
response of the chain in such a manner that the persistence length appears effectively reduced according to the relation 



For a single loop this is obviously a finite size effect involving the scaled total length l p /L, but the effect remains 
significant over a large range of parameters (e.g. l a p pp m 0.7 'Alp for L/l p — 50). Therefore the interpretation of 
corresponding stretching data has to be taken with care: Even though the data seem to suggest WLC behavior, the 
extracted value of persistence length might not reflect the real chain stiffness. 

An intriguing example that actually inspired this work is the force-extension characteristics of DNA in the presence 
of condensing agents like spermidine or CoHex [LH . It shows in some cases a stick-release pattern which might be 
attributed to the sequential unpeeling of single turns of a toroidal condensate [2f| . What is important here is that 
in between the force peaks one can nicely fit WLC behavior but the persistence length that one extracts from these 
data is typically much lower than that of DNA. Only when the last turn is disrupted and the DNA is in a straight 
configuration one finds the expected value of the chain stiffness. 

Before going into any kind of theoretical analysis of looped DNA under tension it is important to note that such a 
configuration is intrinsically not stable and has therefore to be stabilized by some mechanism. Some possible mecha- 
nisms are listed in Fig.^ ( a ) supercoiling in twisted DNA (the same phenomenon like in a looped telephone cable), 
(b) DNA adsorption on a surface (e.g. a liquid membrane) [2(|, (c) DNA in a dense liquid crystalline environment 
kinetically prohibiting the loop unfolding, (d) DNA in a strong magnetic field that tends to align it in a plane per- 
pendicular to the field lines |22| and (e) DNA condensed by multivalent counterions and other ligands that form a 
freely sliding link. 




(1) 



3 



In Section II of this paper we will shortly review the Euler-Kirchhoff elastic description of the (constrained) ground 
states of DNA under tension. It is extremely useful for understanding the behavior of constrained "cold DNA". 
By "cold DNA" we metaphorically mean DNA in situations where the importance of its configurational entropy is 
negligible as compared to its elastic energy. This is typically the case for short DNA lengths (below its persistence 
length lp) and large energy densities (larger than tens of fcgT's per lp). In the second subsection we switch on 
the temperature and discuss how the thermal DNA wiggling affects its behavior. As mentioned above such "hot 
DNA" responds purely entropically to moderate pulling forces. We review the well known derivation of its mechanical 
"equation of state", i.e., the force extension behavior of stretched DNA. In Section III we derive the statistical 
mechanics for looped DNA under tension for the most simple case where the looped DNA is confined to two dimensions 
as depicted in Fig.^3. In this context we will learn how stretched DNA behaves when its new "ground state" is far 
from the straight configuration. The analytical machinery that is applied and developed further here has its roots in 
classical problems of physics like instantons in quantum mechanical tunnelling p^ . The unifying concept behind all 
these phenomena is that of path integration in the semiclassical limit. In Section IV we finally calculate the stretching 
of looped DNA in three dimensions. We start with the case of DNA being oriented in a strong magnetic field, Fig.^i. 
After having given a rigorous derivation of this case, we determine the partition function for a loop stabilized by a 
sliding ligand, cf. Fig. ^ which finally leads to Eq. ^ 



II. INTRODUCTORY REMARKS 



A. Euler elastica 



The basic assumption of a purely elastic description of DNA (and other semiflexible polymers as well) is that the 
local energy density of a given DNA state is given as a quadratic function of the underlying distortions from the 
straight state. Let us consider the simplest situation where the DNA twist degree of freedom can be neglected. This 
can be done in cases when the DNA twist is not constrained from outside, i.e., when no external torsional torques are 
acting on it. Then we can describe the path of the DNA of given length L and bending constant A subjected to an 
applied tension F by the space curve r (s) with the tangent t (s) = J^r (s). It is convenient to choose the parameter 
< s < L as the contour length and to normalize the tangent to unity t 2 (s) = 1. The elastic energy under an applied 
force F writes in this case 28] 



•«-r(4(s 



P • t \ds 
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Here A is the bending stiffness that is usually expressed as A — Ipk^T where lp is the orientational persistence 
length; for DNA l P « 50 nm [29|. 

Let us look first at "cold" DNA, i.e., at a molecule shorter than lp where we can in principle neglect entropic 
contributions to its behavior. The problem of finding the DNA conformation reduces in this case to the classical 
problem of inextensible elastic beam theory (2^| of finding the energy minimizing state SE/St = which satisfies the 
given constraints. In a concrete computation one would parametrize the unit tangent vector t in spherical coordinates 
t = (cos (fjsmO, sin </>sin#, cos 9) and put the force along the z— axis so that the energy now writes 
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Note that this linear elastic ansatz can be readily extended to the description of twisted DNA states [3(j by introducing 
another degree of freedom, the twisting angle ip(s), in addition to 4>(s) and 9(s). In this case one has to modify 

Eq. 01 by adding the term -|- ^</>cos# + ip\ with B denoting the twist-rigidity constant that is for DNA of the same 

magnitude as the bending constant, B « lOkpT nm p9l |. The reason why we can neglect it in some (but by far not 
all) problems is that if the twist angle -0 is not explicitly constrained (no rotational torque or torsional constraining of 
DNA) ip can always adapt so that the B multiplying term in the integral vanishes (without affecting cf) (s) and 9 (s)). 

Remarkably, as pointed out by Kirchhoff [31| the total energy of deformed DNA (elastic rod) can be mapped onto 
the Lagrangian action of a symmetric spinning top in a gravity field. The angles then 6 (s) ,</> (s) and ip (s) describing 
the local deformations of the rod along the contour length s become the Euler-angles 9 (t) ,<f> (t) and ip (t) of the 
spinning top describing the rotation of the internal coordinates system (with respect to the space fixed frame) as 
functions of time t. All the quantities appearing in Eq. [3] have their counterparts in the spinning top case |32j |. The 
tension F is the equivalent of the gravity force acting on the spinning top; the rigidity constants B and A correspond 
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to the principal moments of inertia around the symmetry axis and perpendicular to it, respectively. The resulting 
rod shapes are usually called Euler-Kirchhoff filaments (in 3D) or Euler-elastica (in the 2D case). They are given 
explicitly in terms of elliptic functions and integrals [23 . The latter fact allows one in many cases (for a given set 
of forces and boundary conditions) to obtain the DNA shapes in an analytical or at least numerically inexpensive 
manner. 



B. Semi-classical straight DNA stretching 



The previous description of DNA conformations via the ground state of a purely elastic beam can, however, only 
be successful for very short DNA (shorter than its persistence length lp). In many practical situations with the DNA 
molecules having lengths on the order of microns to centimeters (3> lp) one needs to go beyond the ground state 
description. But as shown in this paper, this is even true for short DNA (< lp) forming a loop. An important question 
(from the experimental and theoretical point of view) in the context of DNA stretching is the determination of the 
mean end-to-end distance of a DNA chain as a function of the stretching force F at a finite temperature T [lj ■ Then 
in order to take into account temperature effects on a WLC under tension one has to compute the following partition 
function 

Q(F,L, T) = f S (t 2 - 1) V 3 [t] e~P E ® (4) 



with the energy expression given by Eq. [21 This partition function is formally the imaginary time analytical continua- 
tion of a path integral of a quantum particle on a unit sphere subjected to an external force. The chain inextensibility 
constraint (represented by the S function in Eq. H} makes this path-integral a highly non-trivial quantity to compute 
as it introduces a parametrization-dependent nontrivial measure term. But as shown below by choosing a proper 
parametrization of the unit-sphere this unpleasant term does not give any contribution if we limit the computation 
to the semiclassical approximation. 

Let us briefly rederive the well-known results for the force-extension behavior |l2j |. To evaluate the path integral 
we introduce the following representation of the tangent vector defined in a Cartesian coordinate system [2^ i = 
(y 1 — q 2 , q y ,q z )- If we now consider a force that points in the x direction, then the energy becomes only q dependent 
and reads 



with = y, z). The metric and its determinant that determine the 0(3) invariant measure are respectively gij (q) = 
$ij + i—g2 an d .9 (q) — a , and the partition function is now written in a curved space as 

Q(F, L,T) = J V 2 [q] Vg¥)e- f3E[q] = J V 2 [q] e -^M-*»M (6) 

with a measure term that can be exponentiated E m [9] = L dt log (l — q 2 ) which is highly singular. To go further 
we now parameterize the unit vector by the Euler angles (j)(s) and 9(s). The choice of the force direction parallel to 
the ir— axis in Eq. [S] turns out to be necessary because an expansion of Eq. [3] around the straight configuration 9 = 
is singular as the angle <fr(s) is then arbitrary. This causes no technical problem when dealing with ground states 
since <f> enters Eq.[3]only through its derivative <p but we need to rotate the force direction into the a;— direction before 
dealing with the statistical mechanics of "hot DNA" on basis of Eq. Introducing the angle #(s) = 9(s) — tt/2 one 
has t = (cos cos i9, sin i)(> cost?, sini9). Having J = cos 2 i9cos^ as the Jacobian of the transformation the partition 
function writes: 



q(f,l, t) = I v[<j)]v [§] e -mm-E m m ^ 



with the elastic energy 



L / A f 



E [#,(/>]= J 1^ ( </> 2 cos 2 ?? + i? 2 ) +Fcos?9cos0) ds (8) 



and with a measure expressed as 



E m [4] = -6(p) [ dslog(|costf|) 
Jo 



(9) 
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which guarantees the 0(3) invariance of the measure V [(f)] T> [d] of the functional integral. Our "non-standard" 
parametrization of the unit vector tangential to the chain and the choice of a force pointing in the x-axis that 
looks unusual are necessary in order to deal properly with the measure in a semiclassical approach of the nontrivial 
functional Eq. [7| Instead, the standard trick for WLC is based on an analogy between the partition function and the 
Feynman amplitude of a quantum particle. The partition function is then approximatively evaluated by determining 
the eigenstates of the associated quantum Hamiltonian [T2I ITsI ] . This method seems to be difficult to adapt in the 
presence of non- trivial saddle points, even though it has been applied for tightly bend DNA [33 • 

Computationally Eq.[S]with the two functions t9 and (f> entering the energy in a nonlinear manner makes the problem 
difficult to be treated analytically. We therefore use the harmonic approximation valid for small fluctuations around 
the straight configuration (t || e x ), i.e., we expand the energy[S]at the quadratic order around the trivial saddle point 
(<f>o = $o = 0). As (3 is absent in front of the measure, Eq. this later does not participate to the selection of the 

saddle point, but one has to take it into account when considering quadratic fluctuations, i.e., E m ps —5 (0) J Q L dsd^. 
The saddle point being trivial the measure term vanishes and the partition function factorizes into two independent 
partition functions: 

Q(F,L,T) = e^ FL Qj(F,L,T) (10) 

with 

Q 1 (F,L,T) = j ' V[^]e-^o{M 2 +F^)ds (n) 

Note that this factorization property and, in particular, the cancellation of the measure are due to our choice of the 
coordinate system. 

In order to compare later the free energy of the straight chain with that of the looped configuration, we compute 
the path integral with the boundary conditions <p(Q) = d(0) = and 4>(L) = fl(L) = which are the most convenient 
choice for a semi-classical evaluation of the path integral around a non-trivial saddle point. The Fourier decomposition 
is then restricted to sine functions 4>{s) = ^2 \/2/L sin (co m s) 4> m with frequencies uj m = irm/L. The evaluation of 
path integral then reduces to the computation of a product of Gaussian integrals leading to 

R^fTF e /3FL 



Q(F,L,T) = ^- — (12) 

27r sinh (lM 



The force-extension relation can then be deduced from the expression (Ax) = —dG/dF where G(F,L,T) = 
— 4lnQ(F, L, T) is the free energy of the system. We then obtain 

A '" ' -^=coth(L/A) (13) 



L 2/3FL 



Here we introduced the quantity A = yj A/F, usually called the deflection length or tension length \2(X l34| , that 
becomes the relevant length scale in the case of DNA under tension replacing the usual (tension-free) persistence 
length l P = A/k B T. 

^From the force-extension relation Eq. 1131 we see that two limiting cases corresponding to regimes of small forces 
L/A<1 and large forces L/A> 1 can be studied analytically. 

1. Small forces regime: L/X <C 1 

One can readily see that this condition implies a small force regime (3FL <C l p /L which is compatible with the 
harmonic approximation only if the persistence length is much larger than the chain length (l p 3> L ). Then in this 
case Eq. ED becomes 

(Ax) L L 3 F L 

^ '- fa 1 + ps 1 (14) 



6l p 24lp A 6l p 



(15) 



i.e. thermal fluctuations lead in leading order to a force-independent small reduction of the end-to-end distance. 
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2. Large forces regime L/X 1 



This regime implies the condition (3FL 3> l P /L, that can be made compatible with the harmonic approximation 
for any value of the ratio l p /L. The free energy of the WLC under tension is then approximately given by 



G(F,L,T) 



2tt 



whereas the force-extension relation in this limit becomes 

(Ax) . I 



1 



L 2PVFA 2f3FL 

In this force regime the term 0(1/ (3FL) can always be neglected and one arrives at the important formula |l2j | 

(Ax) k B T 



(16) 



(17) 



(18) 



This can be solved for the force: 



F 



4Z, 



(l-(Ax)/LY 



(19) 



The force, Eq. 1191 is of entropic origin as the proportionality to temperature indicates. Equation 1191 turned out 
to be a very powerful tool for directly and accurately determining the persistence length of DNA molecules from 
micromanipulation experiments under a multitude of conditions [2j • One should note that Eq. 1191 is only valid in the 
limit of large forces (F ^> ^f^- = 20 fN) and large relative extensions (Az) /L « O(l). Looking at its simplicity it is 
somehow surprising that it is experimentally accurate for piconewton forces almost up to the point where DNA starts 
to melt and the WLC description breaks down (around 60 pN). 

To have an expression that also includes very low forces (on the femtonewton scale) one usually uses the following 
interpolation formula |12| 



F 



4 1 



L 



(Az) 
L 



(20) 



In the limit of small extensions, (Az) /L <C 1, one recovers F — | kB J^ z ^ which is the force that one expects for a 
Gaussian coil perturbed by weak forces . For large forces one asymptotically recovers Eq. 



III. THE LOOP IN 2D 



In the following we consider a DNA chain under tension that contains a sliding loop. The corresponding shape 
at zero temperature is that of the homoclinic loop 361 that is a solution of the Euler-Lagrange equations. This 
filament shape that was already considered by Euler [32j | is two-dimensional. For any given finite tension F the 
homoclinic loop turns out to be stable for arbitrarily large in-plane perturbations. Indeed the 2D homoclinic loop 
can be considered as a (static) topological soliton appearing in many contexts of contemporary physics ranging from 
Josephson-junctions, dislocations in solids to QM tunneling problems j^. In the current section we study the 
DNA chain being confined to two dimensions; only in Section IV we go into the third dimension by allowing also 
out-of-plane fluctuations. The problem is then that the loop is intrinsically unstable (in contrast to a false claim in 
literature 38j) and one has to introduce potentials or constraints necessary for its stabilization. 



A. The partition function 

Consider a looped DNA chain under tension F along the x-axis. In this section the DNA is only allowed to 
fluctuate in-plane (as it is the case for a chain adsorbed on a fluid membrane, cf. Fig. 2J)). We neglect the DNA 
twist degree of freedom that in general - if not explicitly constrained - decouples from the DNA bending energy. 
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To obtain the force-extension behavior of the loop in 2D we evaluate semiclassically the partition function Qi oop by 
considering quadratic fluctuations around the saddle point 4>i oov that is here the loop configuration. We impose that 
the angles at the extremities of DNA are clamped in an orientation parallel to the force direction, so that <j>{—L/2) = 
and <p(L/2) — 2ix. Then the partition function in 2D corresponds to the following quantum probability amplitude 
expressed in terms of a path integral: 

j.(2n,L/2\) 

Qioop = (0, -L/2A|2tt, L/2X) = V[4>] e~ f3E ^ (21) 

J(0,-L/2X) 

where the energy can be written: 

, r L/2X (\ . \ 

E [4>] = VAF / \-4> 2 - cos cj> dt (22) 
J-L/2X V 2 / 

with the dimensionless contour length t — s/A; dots represent from now on derivatives with respect to t. In the spirit 
of the Kirchhoff kinetic analogy from Section IIA the bending energy in Eq. [221 corresponds to the Lagrangian of 
a spherical pendulum in the gravitational field. We now expand E [<j>\ up to quadratic order around the minimum 
configuration 4>i oop by introducing a fluctuating field 8(j> such that <j) = 4>i oop + d(j). The linear term SE in this expansion 
vanishes because 4>i 00 p is an extremum point of E and we have 

E [()>ioop + 5(j>] = Ei 00 p + 5 2 E [4>ioo P ] (23) 



1. The saddle point contribution 

To determine the saddle point configuration we solve the Euler-Lagrange equations of Eq.[221that gives the following 
nonlinear equation 

4> = sm<i> (24) 

which is the time independent Sine-Gordon equation well known and studied in many systems especially in the context 
of solitons (and their applications like Josephson junctions, cf. Davydov's book |23)- Beside the trivial solution cf> = 
that corresponds to the ground state but cannot describe a loop configuration there exist other topological solutions 
of Eq. [2] that are appropriately called solitons or kinks. Now Eq. [21 can be integrated twice to obtain 



J<t>{t ) V 2 [C - cos <j) ') 

with an integration constant C. The general solution of Eq. EH with arbitrary C leads to elliptic functions. With the 
condition to = and <f> (0) = 7r the solution reads 

,2/ ' 



m 



cos <j)i 00 p (t) = 2sn 2 _| m - 1 (26) 



(f>ioo P (t)=%+ 2am I — =|m ) (27) 



m 

with sn and am being the Jacobian elliptic function with parameter m 39] whose value is related to C in Eq. 1251 via 
in = 2 J (1 + C). The parameter m with the range ^ m $C 1 results from the clamped boundary conditions and is 
implicitly given by 

L L 



mK(m) = — = -^F/A (28) 

with K (m) denoting the complete elliptic integral of the first kind |H3 • I n the Kirchhoff analogy the solution Eq. [2H 
describes a revolving pendulum that makes one full turn during the "time period" L/X. 

The " classical" (T — 0) bending energy of the loop as an implicit function of the force is then given by 

PE [4>ioo P ] = 4 l j-K (m) [K(m) (m - 2) + 4B(m)] (29) 

where E (m) is the complete elliptic integral of the second kind |39j . We compute now the contribution of the quadratic 
fluctuations around this looped saddle point to the partition function. 
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2. The fluctuation contributions 

In the semi-classical approximation the partition function Eq. 1211 can be written as a product of an energetic 
contribution and a quadratic path integral over the fluctuating fields 8<p satisfying the Dirichlet boundary conditions 
ty(-£) = ty(-&)=0: 

Q loop = e -^ loopQ fluct m 

with the partition function corresponding to the quadratic fluctuation contributions given by 



>"-' / VMe-^5-^ s ^ = J ft^g , (31) 



2nD(-±± 



2sn 2 [-=\m -1 (32) 



Here D(—j^, j^) is the determinant associated to the quadratic fluctuation operator T that reads 

. _ , / « 

The problem of finding the eigenvalues of this operator falls into a class of "quasi exactly solvable" problems and 
typically appears in quantum mechanical problems. The corresponding differential equation is called the Lame 
equation |40|. It admits simple solutions in terms of polynomials of elliptic functions sn, cn and dn. Its discrete 
spectrum is known |4Cj and writes 

V-\ = 1 and /_i (t) = cn [ |m 

to V Jm 



isq = and Jq (£) = dn ( — ^=|to J (33) 



V\ = — and fi (t) = sn [ \m 

to \ \/ro 



One sees immediately that the only eigenfunction satisfying the Dirichlet condition is f_i so that the smallest eigen- 
value of T is v-\ that we denote in the following by 

1 — m , 

Mo = (34) 

to 

Therefore for a molecule of finite length there is no zero mode but [1q goes to zero in the limit of infinite length L 
that in terms of m corresponds to the limit m — ► 1. The existence of a vanishing eigenvalue is the consequence of the 
translational invariance t — > t + t of the loop solution that formally causes a divergence of Eq. E2 Hi . 

The determinant D (— Jr, can be computed directly via the method of Gelfand and Yaglom |41j that consists 
of solving an initial value problem on the interval [— L/2X, L/2X]. The explicit solution for D (^, —j^) can be stated 
in terms of the "classical" solution 4>i oop (t); 

I' I) = +*» (I) ^ (-<§) j_ L/2X jrf—^ 

= 2^iE (to) (35) 
This expression, however, has to be taken with caution since it is incorrect for large values of L/X missing a factor 

L / F 

~ e~~V~x that corresponds to the fluctuation contribution of the linear part of the DNA. To solve this problem one 
has to take the translational invariance of the loop into account. The way to deal rigorously with a zero mode in the 
infinite L case is well known [22^. One has to consider the infinite number of degenerate saddle points resulting from 
the translational invariance by considering the collective time coordinate to as an integration variable instead of the 
normal mode do associated to the zero mode eigenfunction fo(t). For a finite chain length we also have to take into 
account the displacement of the kink solution but this time only on a finite interval of length L. We compute the 
corrected determinant D corr by removing the would-be-zero mode from the determinant and by considering explicitly 
the finiteness of the space that this mode can populate: 



D - [ -TyTx)) = ne -^ ^ ao { D {-2-x>2x)) (36) 
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where the boundary a must be determined by computing the Jacobian defined by Sciq — J 1 (m) £to- To do this 
consider a small translation of the loop which is then given by 8<pi 00 p — i'loop&to ■ Because the same translation can be 
done by the cigcnfunction fo(t) we have 5<fii oop — fo(t)Saa. It is easy to check that the normalized-to-one eigenfunction 

/ L /2A * \ — • 

associated to the zero mode is of the form f (t) = I j_ L / 2 \ dt4>j oop ) <filoop(t) from which follows that the Jacobian 
is simply equal to the normalized factor 

■'-»=(/ =(^&-Y 2 (37) 




^From this relation we deduce a — ^ ( ^' E j^ l ^j and finally obtain the partition function 



Q \ erf J ^pm k e -^mmKm- Wm)] {m) 

\ AiTy/mE (m) / 1 | %Jm A / 

Using the relations 1281 and l3~4l we can rewrite this expression fully in terms of m: 

Qioo P (m) = fefg) 1 erf ( 2 ( 2±) " if* (m) E* (m) (1 - m) A e -4£*M[*( ro ) (m -2) + 4 EM] (39) 



Note that the erf-function only differs significantly from unity for values ofmssl which corresponds to the long DNA 
limit L/X 3> 1, i.e., the correction given by the entropic contribution of the loop is significant only in this limit. This 
is important because our computation of the Jacobian is strictly valid only is this limiting case: The eigenfunction 
fo(t) associated with the zero eigenvalue v$ (and responsible for the translation of the loop) satisfies the Dirichlet 
condition only in the limit of an infinitely long chain. This is why the zero mode (and hence fo{t)) are excluded from 
the determinant for finite chains. Physically the boundary condition used in our computation implies that for a finite 
chain the shift of the loop induces an elastic deformation that costs energy. Only in the long-chain-limit the loop 
can move freely. The constraint Eq. [2H|can be solved for m w 1 giving m 1 — 16e~~ showing that the smallest 
eigenvalue ^iq, Eg. 1341 reads /iq = 16e~~ . So indeed - as intuitively expected - this eigenvalue becomes asymptotically 
zero for L/X — > oo (i.e. m — ► 1). 



B. The force-extension relation 



The force-extension relation of the looped chain in 2D follows from the free energy Gi oop — — In (Qi oop ) via 
(Ax) = —dGioop/dF with Qi oop given by Eq. 1391 Due to the structure of Qi 00 p the mean extension is a sum of three 
terms: a contribution from the bending energy, a second one from fluctuations around the loop configuration and a 
third from the error function, i.e. 

(Ax) = (A E x) + (A||x) + (A err z) (40) 
The saddle point contribution to the mean extension is given by 

(A E x ) = — —— 41) 

m K (m) 

The contribution resulting from fluctuations around the loop configuration (determinant) is given by: 

x _ L 2 E 2 (m) + (1 - m )K*(m)- 2(1 -m)E(m)K(m) 
\ 11 ; m p m [ (E(m)K(m)) 2 1 [ ' 

Finally, the contribution coming from the error function is: 

L(l- m)[3E 2 (m) - (1 - m)K 2 (m) - 2E(m)K(m)\ e -^(i- m )E( m )K 3 ( m ) 
{i\ crr x) — ■ - J — . 

2V2^merf ^2 (2^) 2 K% (m) E^ (m) (1 - m^j E(m)<J l f(l - m)E(m)K(m) 

This allows one to immediately plot force-extension curves for a loop in 2D. We dispense here with giving such a plot 
since the curves turn out to be very close to the corresponding ones in 3D, presented below in Fig.^J Instead we only 
extract from Eqs. ^2 to 021 the limiting cases of small and large forces. 
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1. Limit of small values of m w 0. 

This corresponds to a regime of small forces L/2X <C 1, valid only for chains satisfying l p /L 1. In this limit the 
functions sn(x|m) ~ sin(a;) and am(x|m) ~ x, so that the loop configuration given by Eg . 1261 corresponds to a circle 
(f>ioop (s) = n + 2stt I L + 0{m). The bending energy is then given by: 

pE chcle = 2tt 2 j- + 0(m) (44) 

where we used K(m), E(m) « n/2 + 0(m) for m very small. 

To determine the force-extension relation we expand the various contributions in expression Ea. l41l to the first order 
in m and replace m by =y £ . We then arrive at 

< a ->^i^ + °<" i2 » < 45 ' 

which shows that the bending energy contribution to the elongation goes to zero with the force. This is expected as 
the bending energy is independent of the force when m goes to zero by virtue of Eq. 1441 

In the same manner we obtain for the contribution due to the quadratic fluctuations around the loop (there is no 
translational invariance in this case as the linear part of the chain is negligible) the expression 

(A F )« i ^-(l + 0(m)) (46) 

We find here that the thermal fluctuations cause on average an increase of the end-to-end distance (resulting in a 
reduction of the loop size) . Note that this is contrary to the stretching of a linear DNA where entropic effects lead to 
a shortening of the polymer (cf . Eq. I14|) . 

As the loop cannot slide in this context, the contribution of the error function to the force-extension relation should 
vanish. Indeed our computation gives 



120 L 2 

which is negligible due to l p /L ^> 1 



(A err x) » - — — e~ 2L (l + 0(m 2 )) (47) 



7T 2 l p 



In conclusion in the regime |- <C 1, the force extension relation is given by: 



{AX) ~ L -(l + ^) (48) 



L 4tt% v 2 



i.e., for short looped chains the extension grows linearly with the force. 



2. Limit m « 1 : the homoclinic loop case 

In the limit m — ► 1, K(m) diverges as In (4/yI — m) and E(m) m 1. By virtue of Eq. 1281 this corresponds then 
to the case L/2A^> 1 (strong force regime) where the length of the molecule is very large compared to the loop size 
of order A. If one is only interested in the force extension curve one can directly take the limit m — > 1 in Eq. 1401 
It is, however, interesting to rederive it via the saddle point approximation of the path integral in the infinite chain 
limit because then the strong analogy between our computation and the semi-classical treatment of the tunneling of 
a quantum particle in a double well potential becomes very transparent. 

In the limit of a very long DNA chain, expressions l28l reduces to a kink configuration 4>i 00 p — 4 arctan e* interpolating 
between the two values <j>{—oo) — and </)(+oo) = 2-7T, cf. also Fig.|2 Equation 1261 is then given by 

2 

cos (j)i 00 p (i) = 1 -s— ■ . (49) 

cosh (t) 

This saddle point solution is correct only in the infinite chain limit, but for finite large length the corrections are of 
order e~ L / A . This implies that the bending energy of the kink is then given by 

Eioop = E [4> loop \ = -FL + 8VAF + O (e- L/x ) . (50) 
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The fluctuating quadratic operator in this long chain limit is 

d 2 



ofi ' v'-JrfwJJ (51) 

This operator is the same as the fluctuating operator obtained by considering fluctuations around the kink solution 
that appears in the semi-classical treatment of a quantum particle in a double well potential. Having this operator 
one can compute its set of eigenvalues by solving the Schrodinger equation for a particle moving in a potential of the 
Rosen- Morse type (see |2^,|4^). Then the partition function is given by 



Q loop » i/3LFe-*vT e H»(8^-"0 (52) 

from which we deduce the free energy 




Gioop -^]Jj+ 8 VAF - LF - i In (^PLF) (53) 

We compare now this free energy to that of the straight state. Note that we cannot use Eq. QUI since it corresponds 
to 3D case. Instead we need the 2D free energy that derives from a partition function that is evidently given by 
Q (F, L,T) = e l3FL Qi (F,L,T). This means that the second and third term of the free energy expression, Ea. 1161 
have to be divided by 2. Subtracting that free energy, Go, from Gi oop leads to the free energy difference: 

AG ioop _ = G loop - Go = 8VAF + o(\ In ( (54) 



7 loop— ^loop 

Then we see that the free energy difference AG; oop _o is dominated by the elastic energy part 8V AF which is the 
second term in Ei oop , Eq. 021 The first term —FL is already present in the straight DNA case and cancels in the 
difference. Besides that (typically very large) term there is merely a logarithmic correction. We note that a weak 
coupling of the thermal fluctuations to a DNA shape has also been observed by Odijk |3^| for circular DNA rings. 
The force-extension curve of a 2D loop is then calculated via (Ax) = —dGi oop /dF: 

^-1-^+4^ (55) 



L V 4 LJ V^FTp f3FL 

To understand the origin of each of the various contribution to Eq. 1551 we consider now the limit m — > 1 of Eq. 1401 
We find for the energetic contribution 

^K^^-iS) (56) 

The contribution of the fluctuation around the saddle point (the determinant) is given by 

(Aiix) w — !— (57) 
x 11 1 m-l A(3F v ' 

We have already seen in the case of the linear DNA stretching that such a contribution is negligible in this force 
regime, cf. Eq.^l That means that in this regime the contribution of the fluctuations to the force extension relation 
around the loop are negligible in comparison to the contribution coming from the elastic energy of the loop. 

Now the error function must play an important role as the entropy of the loop is no more negligible for a very long 
chain. The contribution from the error function gives 

T Q 

(58) 

The first term in this equation corresponds to the fluctuation of the linear part (in 2D) of the DNA (compare with 
the corresponding 3D-term in Eq. I17f) whereas the second term is negligible. Combining the different contributions 
we recover equation Eq. 023 We may drop the last contribution ksT/FL that is for all practical purposes negligible. 

We can now compare the "equation of state" of the looped DNA, Eq. [23 with the one for the straight configuration 
in 2D given by {Ax ) = -dG /dF: 
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FIG. 2: a) The definition of the Euler angle cj> and the scale of the loop. The loop head diameter (red) is approximately given 
by A = \J A/F b) The loop solution 0; oop (s) as given by Ea. 1491 



Comparing Eq. 1551 and Eq. 023 we see that both have a leading term proportional to F -1 / 2 ; only the prefactor in 
Eq. 1551 is renormalized by a contribution stemming from the elastic part of the loop free energy. 

This implies a fairly simple prediction that is useful for the interpretation of experimental data: Suppose one 
performs a single molecule stretching experiment with a DNA chain that contains a loop. If one does not know about 
the presence of the loop one will fit the data by the usual WLC expression, Eq. 023 and is happy that it works well 
(at least up to the leading term F -1 / 2 ). From that fit the total length of the DNA is recovered correctly (from the 
asymptotic line on the Ax axis) but something strange seems to have happened to the "persistence length" - it is 
smaller than expected. The explanation is simple: The apparent persistence length becomes 

IT = , ~"TT2 ^ 2D (60) 
(1 + 16$) 

This formula is similar to Eq.Q] announced in the introduction, with a factor 16 instead of 8 in front of the Zp/L-term. 
The difference comes from the fact that we allow here only fluctuations in 2D. The 3D case will be studied in chapter 
IV. 



IV. THE HOMOCLINIC LOOP IN 3D 



After having understood the behavior of the homoclinic loop in 2D it seems that a generalization to the third 
dimension should be straightforward. But as we will see there are several traps and some interesting physics on the 
way. The first and main problem is the fact that the homoclinic loop is (unlike in the 2D case) elastically unstable. A 
simple way to see this is to take an elastic cable, make a loop in it and to pull on it (without torsionally constraining 
the ends). Only if we force the loop to stay in a plane (for instance, its own weight can perform this task if the cable is 
lying on a table provided that we do not pull too strongly) it represents a topological excitation that cannot leave the 
rod (except at either of its two ends). So if there is any interesting physics of 3D homoclinic loops it will have to come 
through constraints or loop stabilizing potentials. In the following we mainly consider two stabilizing procedures: In 
one case we remove carefully the unstable mode from the partition function (the loop is then approximately forced to 
stay in a plane) and in the second case we evaluate the partition function in the presence of an explicit self-interaction 
that stabilizes the loop. It is then shown that for very long chains the fashion via which the loop is stabilized is 
irrelevant with regard to the determination of the force-extension relation. 
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A. The unstable and rotational zero mode 



We first discuss here the relevance of the right parametrization of the unit tangent vector. Beside the fact that our 
parametrization allows us to deal properly with the measure term, the importance of this choice appears even more 
evident when considering the 3D loop. Suppose that we study the equilibrium property of DNA pulled by a force in 
the z-direction. The bending energy is 



E 



[ L/2 \- ( 

J -L/2 U V 



6 2 sin 2 9 + 9 2 ) - Fcos9 



ds 



(61) 



and the saddle point is now 4>i oop = and 0i oop given by Eq. 1491 (with 4>i oop replaced by 0i oop ). Looking at small 
variations 89 and 6<p around the homoclinic loop solution we find a positive definite second variation of the energy 
functional, 



5 2 E = 



L/2 
-L/2 



^S9 2 + ^ cos (9 loop ) 69' 



Asin 2 {Oioop) 6<j> 2 



ds 



(62) 



This is in striking contradiction to the expected elastic instability of the loop in 3D. The reason is that the coordinate 
system has a singularity at = where the </>-angle becomes arbitrary. As explained in section IIB, the way to 
circumvent the problem is to rotate the force direction and to put it along the x-axis so that the potential energy 
part writes now — F cost/) sin 9. In terms of the angles <f> and *& = 9 — tt/2 the energy writes now 



E [■&, (f\ = VAF 



I-L/2X 



J-L/2X 

with the corresponding Euler-Lagrange equations 



i) 



4> cos ■& — 2(j)d cos d sin i? 



cos 2 i? + v ) - cos 6 cos d dt 



cos 4> sin i? 
sin 6 cos ■& 



(63) 



r cos § sin d 



(64) 



We choose in the following the d = solution, i.e., we put the loop into the x-y-plane. This imposes no restriction 
as we can always rotate the coordinate system around the x-axis to achieve ■& = (). In this case we have <f> = sin</> 
which is the same as Eq. and the saddle point is then given by 4>i 00 p (t) , cf. Eq. |57| and $i 00 p (i) = 0. Then by 
considering small fluctuations around this saddle point that satisfy Dirichlet boundary we can again expand Eq. 1631 
up to second order and obtain 



0E [Si), cj>i oop + , 



(3Ei oop 



[3VAF 
2 



L/2X r-L/2\ 

S(/) (t) f || (t) 5(f) (t) dt+ 5d (t) Tj_ (<) 5i9 (t) dt 

-L/2X J-L/2X 



(65) 



where the loop energy Ei oop and the in-plane fluctuation operator Tii are of course the same as in the 2D case, 
cf. Eq.HDand Eq. 01 respectively. New is in Eq. [HS] the out-of-plane fluctuation operator T_l given by 



' dt 2 



6sn^ 



(66) 



Note that with our choice r di oop (t) = the measure term does not contribute at this level of the approximation. The 
main consequence of the quadratic expansion around the saddle point configuration is that the variables i? and <j> 
decouple so that the full partition function Qi oop factorizes into the product of the 2D partition function Q2D (given 
by Eq. I39fl and the partition function Q± accounting for out-of-plane fluctuations: 



Qloop — Q2dQ± 



(67) 



Although very similar to Ty the behavior of the out-of-plane operator Tj^ is fundamentally different. The discrete 
spectrum of Tj_ consists of two eigenvalues fi J z 1 = — 3/m and [1q = 0, the first of which is indeed negative |40| . 

The zero eigenvalue mode of Tj^ comes from the rotational symmetry around the x-axis in a similar manner as 
the translational invariance of the loop causes a vanishing eigenvalue of Ty . To compute the contribution of the 
infinite number of degenerate saddle point related by a rotation around the x -axis we look at infinitesimal rotational 
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transformations of the loop in 3D. It is straightforward to show that up to quadratic order a rotation of a kink with 
$ioop — around the rc-axis by a small angle Se corresponds to the following small changes in "&i oop and 4>i 00 p'- 



Sfiioop ~ Se sin (pi oop 
-Se 2 smd 



™loop 



'loop 



COS( 



Hoop 



= O (Se 2 



We note that in lowest order this rotation leaves 4>ioop unaffected, so that the same rotation S"di 00 p ~ Se sin 4>ioop 
can be done also by the normalized eigenfunction 'doit) associated to the zero mode ($ and <f> formally decouple) 
alone, 8-di oop « Sbo'do (t) where bo is the normal mode variable associated to i?o(i). It is easy to check that this mode 
normalized to one writes 




-1/2 



<>i> (I) = | / <// hii "J,,,,,, } sm4>i oop 
where sin0/ oop = cn ( -|=|to) sn (^-j=\rnj . Then the Jacobian defined by dbo = J~ 1 de is given by 



(68) 



J- 1 



*L/2\ 



1/2 



I-L/2X 



dt{ddi oop /de) 



r L/2X 



1/2 



dt sin 



I-L/2X 



"loop 



or explicitly 



/o\V2 

J- 1 (to) = 2m- 3 / 4 - ((2 - m)E(m) - 2(1 - m)K{m)) 1/2 



(69) 



(70) 



This relation will be necessary for the computation of the out-of-plane determinant. In quantum mechanics the ground 
state wave function has no node, the first excited state wave function has one node, etc. In our case the wave function 
i?o has one node so it cannot be the ground state and a wave function with no node in the interval considered must 
exist. It is obviously the eigenfunction of the unstable mode that is given by the (unnormalized) expression 



i?_l (t) = cn 



t 



dn 



t 



with the eigenvalue 



Mix 



(71) 



(72) 



This negative eigenvalue makes the 3D loop mechanically unstable. An overview of the three discrete eigenmodes is 
provided in Fig. [3] 

In the following we will consider three different mechanisms to stabilize this mode: In the next section we study 
looped DNA in a strong magnetic field that breaks the rotational invariance and then we enforce the cancellation of 
the unstable mode either by a geometrical constraint (Section IVB) or by an explicit self interaction (Section IVC), 
both of which keep the rotational symmetry. 



B. DNA in a strong magnetic field 



A physical situation in which the DNA loop is stabilized is if we switch on a (very strong) magnetic field along the 
z-axis perpendicular to the force direction along the rr-axis, cf. Fig. ^1. The DNA nucleotides (having 7r-electrons) 
are known to prefer alignment perpendicular to the field, i.e., DNA exhibits a negative diamagnetic anisotropy |27| . 
The application of a magnetic field H along the z-axis drives the DNA molecule into a plane parallel to the x-y plane. 
The total energy of the DNA writes in this case 



E[4,4>] = 



fL/2 r A 

/ -(</> 2 cos 2 tf + i? 2 -.Fcos, 

J-L/2 L 2 V 



■ cos d + — sin 2 i9 



ds 



(73) 



The last term gives the coupling between the DNA tangent and the magnetic field H where k = —x a H 2 /h characterizes 
the coupling strength. Here \a denotes the (experimentally accessible) diamagnetic anisotropy of a single DNA 
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FIG. 3: The three discrete eigenmodes in action for m « 1: a) The translational mode <f>o = 1/ (\/2coshf) (Eg. I33B . b) The 
rotational mode i?o = \/3/2 sinhf cosh -2 i (Eg. I68H and c) The unstable (out of plane tilting) mode #_i = \/3/4 cosh -2 t 
(Eq.EU- 



basepair (23 and /i = 0.34 nm is the distance between the subsequent DNA basepairs. Note that \a is negative here, 
i.e., k > so $ = is the preferred rod orientation for large K. 

Expanding E [d, <f>] up to second order we obtain the same expression as Eq. 1651 except that Tj_ is replaced by a 
new out-of-plane fluctuation operator T"^ 



The spectrum of Eq. EH is given by shifting the spectrum of Tj_ by the constant k/F leading to the eigenvalues 
/xj = jitj- + k/F. The rotational mode is immediately destroyed for any non-zero coupling constant k > 0. More 
importantly the previously unstable mode $_i now becomes stable provided that k/F > 3, i.e., for k > K cr n = 3F. 
As the partition function factorizes into the product of the 2D partition function and an out-of-plane contribution 
we just need to compute the determinant associated with the fluctuation operator 1741 which is given in Appendix B, 
Eq. IB12| for the case L/2X ^> 1. From this we obtain 



Q ^-V tt( C -2)( C -1) 6 (?5) 



with c = vf + k/F. The free energy /3G = — In (Q2dQ±) with Q2D given by Eq. [52] has the following form 

/3G = —f3FL + WVFA + (1 + v/TT ^ ) - In ( j^g^^ F^ (76) 

Differentiating this expression with respect to F leads to the force-extension relation for all forces F < k/3. Since 
this turns out to be a lengthy expression we give here only the result for the limiting case k» F (and - as assumed 
above - L/2A > 1): 




k B T ( 1 l P 1 /F\ 3/2 „/F\ 5/2 \ „/ 1 



^ V 4 +4 f-8UJ +0 UJ ) + °{Wl) (7?) 

This expression is similar to the 2D-case, Eq. which is related to the fact that we assume a strong cost for 
out-of-plane fluctuations by setting k 3> F. The major contribution from the out-of-plane fluctuations is the second 
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term on the rhs of Eg . 1771 that describes an effective F-independent shortening of the contour length. The next-order 
K-dependent correction is already by a factor F/k smaller and therefore negligible. 

Finally, could we experimentally observe the force extension curve derived above? Unfortunately the coupling 
parameter k turns out to be too small for reasonable magnetic fields |43l | to be of physical relevance in practice, i.e. 
k <C 3F. Nevertheless, the formal diamagnetic term Ksin 2 $ introduced in Ea. l73l is conceptually useful to understand 
the (otherwise unstable) behavior of the DNA loop in 3D. It also turns out to be technically convenient to use an 
infinitesimal small " diamagnetic term" in order to break the rotational symmetry of the system for the computation 
of the 3D determinant when dealing with the rotational zero mode (cf. Appendix A). 



C. Force-extension with a geometrical constraint 



In this section we compute the partition function of looped DNA by forcing the mean tangent of the loop to stay 
in a plane which is the simplest stabilizing procedure. This geometrical constraint corresponds to applying forces at 
the two chain termini that maintain them in-plane. The constraint is implemented by the introduction of a delta 
Dirac distribution in the partition function, so that the out-of-plane partition function in the presence of a external 
magnetic field becomes 

Ql (0 C ) = / S I ~ / Stidt I e ~ S ~^ E ' I - / ^ s&tlMdt V [5#\ (78) 




The formal presence of the external magnetic field is necessary because the rotational mode $o corresponds to a zero 
eigenvalue and causes a divergence of the partition function Q^ =0 . The problem results from the fact that a rotation 
of the kink around the x-axis costs no energy and consequently the entropic contribution of this state space direction 
seems to diverge (within the Gaussian approximation implied by the saddle point approximation used here). To 
circumvent this problem we employ the following trick. Instead of T_l we use from Eg. 1741 and after performing 
all other calculations we let k — > (note that \ k —q = Tj_). Physically this procedure corresponds to infinitesimally 
breaking the rotational symmetry (around the force direction) and restoring it afterwards in a controlled manner in 
the limit n — > 0. 

Physically it is clear that the main contribution of the mean value of the tilting angle defined by 



A f L/2X 

(cW) = - / 5<&dt (79) 

L J-L/2X 



comes from the unstable mode that induces the large out of plane deviation. The contribution from the rest of the 
eigenmodes is small and stable, so that we can make the following approximation 

(W)«a_i<tf_i)=a_iJJ m (80) 

V 2 Kim) ((i + m )E{m) - (1 - m)K (m)) 1/2 

In this way, we approximate the constraint in Eq. 1781 by a constraint that fixes the mean out of plane deviation 
induced by the unstable mode alone to zero. It means that we relax a bit the constraint in Ea. 1781 bv allowing the 
other modes to induce non-zero mean value of the tilting angle. This contribution however will be small and limited 
by the positive spring constants of the stable out of plane modes. 
It is then straightforward to rewrite Eq. EBI as follows 





Ql = »V« It/ ^T- ^ J x ' *-"-*S(a-i)da-i j ^ 2 ^,_ A A) (8D 

where - with k being small - the determinant is now imaginary. Note that by removing the unstable mode from the 
determinant we have taken into account that the Gaussian integral of the unstable mode is given by 

c-n-ig'/a <*" = 1 



and not by as a naive analytical continuation would suggest |22l l44|. Expression Eq. 1811 shows that one 

can not simply remove the unstable mode from the determinant but one has to replace it carefully by introducing 
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a correct constraint expression in the partition function. For instance, cancelling simply the unstable mode would 
introduce non physical divergences in the limit of very small forces. 

We are only interested here in the limit of zero magnetic field. In order to restore the rotational invariance we have 
to deal with the rotational zero mode by dividing out the would-be-zero mode and replacing it by the real physical 
space it populates. Therefore we have to compute the K-independent partition function Q± defined by 

Q ± k = q VM^) (j- 1 (m) J J v^P & ) Ql ( 82 ) 

where the Jacobian is given by Eq. 

The out-of-plane determinant can be deduced from the Gclfand-Yaglom method which specifies that the determinant 
can be obtained from the solution of the following generalized second order Lame equation: 



(f ± + ±) y(t) = (83) 



The determinant is then given by the relation D r j_(—L/2X, L/2X) = y(L/2X) valid when the conditions y(— L/2X) = 
and y(—L/2X) — 1 are satisfied 41]. A detailed determination of the solution of Eq. [S3|is provided in Appendix A. 
In the small K limit the out-of-plane determinant admits the following expansion (cf. Eq. IA36|) 



Dl = -j ^l_ m) ((2 - m)E{m) - 2(1 - m)K{m)) (84) 

which is negative as it should be because of the presence of the unstable mode. Combining the different contributions 
in Eq. |S21we arrive at the following expression for the out-of-plane partition function: 

2V2 fl p V /2 „, wa n _ , 1/2 f(l + m )E( m )-(l-m)K( m )\ 1/2 



Q± = ^{i) *("0 5/a (i-"0 1 ( - ->^""> ^ "v-v-v j (85) 

The complete semi-classical partition function of the looped DNA in 3D is then given by Qi 00 p — Q2dQ± with Q2D 
given by Eq. |2U 

The mean end-to-end distance has now four contributions 

(Ax) = (A E x) + (A l{ x) + (A CII x) + (A ± x) (86) 

where the first three expressions are given by Eq. EJ] to 031 an d the out-of-plane contribution obeys 

L 2 ( 5(1+ m)E 2 {m)- 2(6 -3m- m 2 )K(m)E(m) + (7 -12m + 5m 2 )K 2 \ 

' ±X ' ~ m p mK 2 {m)E(m) \ (1 + m)E(m) - (1 - m)K{m) J 

With the complete analytical expression at hand it is straightforward to compute force-extension curves, some 
examples for different ratios L/l p can be found in Fig.0] The curves show clearly different scaling behavior for low 
and strong forces. In the small m limit (L/X « 1) - corresponding to the case L < l p - Eq. |H3has the following 
expansion 



«-^-(l + 0(m)) 

167T Z i p 



that has the same scaling as the in-plane fluctuation contribution Eq. 1461 

In the limit L/X — > 00 or m — > 1 the out-of-plane partition function, Eq. 1851 takes the following form 

which leads to the force-extension curve 



(/3VFJ iVVWS (89) 



<A ^"-4V^ + ^ (90) 
When adding this result to the contributions stemming from the 2D computation we finally obtain 



/a \ t 1 kBT t a kBTlp 1 9kBT fm\ 
{AX) =L -2\If^ L - 4 V + ~^F~ (91) 

In the following section we compute the force-extension curve in the strong force regime with a loop stabilized by 
a self- attractive potential. This computation allows us to check explicitly that the force-extension relation is fairly 
independent of the details of the stabilizing procedure and a much more physical loop-stabilization leads again to 
Eq. E2 
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FIG. 4: Force-extension curve of a DNA chain under tension with a sliding loop (thick lines, Eg. I8ri|l and without a loop (thin 
lines, Ea. Il.'-ty for different ratios l p /L as denoted by the numbers close to the curves. Specifically we plot here (Fp/F) 1 ^ 2 with 
F p = ksT/lp against Ax/L. In this representation the curves for loop-free chains collapse in the limit of large forces, cf. Eg. 1181 

D. The DNA Self- Attraction and the Homoclinic Loop 

Now we treat an experimentally relevant case in which a loop is stabilized in 3D: self-attracting DNA. DNA is known 
to effectively attract itself in many solvents despite its strong negative bare charge. Typical situations inducing DNA 
self- attraction are poor solvents (like alcohol, small neutral polymers like PEG), the presence of multivalent countcrions 
(like CoHex and Spermidine) or small cationic proteins acting as linkers between two DNA surfaces. Indeed it was a 
single molecule stretching experiment on DNA condensed with multivalent counterions that made us think about 
the force response of loops. 

How should we deal with the DNA self-interaction? A formal treatment that first comes to ones mind is to introduce 
a potential V(||x(si) — x(s2)||) acting between any pair of points si and S2 on the DNA molecule and to write the 
total interaction energy in form of a double integral (over s\ and Sa) as an additional term in our Hamiltonian. The 
problem is, however, that we describe the DNA conformation here by the two spherical angles ($ and cj>) of its tangent 
vector whereas the self- interaction acts in real space ("integrated tangent space"). This makes such a Hamiltonian 
virtually intractable and hence we need a reasonable simplification of the DNA self-attraction. 

To this end we make here two simplifying assumptions: (i) There is only a single discrete DNA self-contact point, 
given by the crossing point of the homoclinic loop solution, (ii) The interaction potential V(||x(si) — x(s2)||) is 
short-ranged enough so that the interaction energy at the crossing becomes independent of the crossing angle, i.e., 
other parts of the DNA (apart from the crossing point) do not interact with each other. 

These fairly reasonable assumptions imply that the loop ground state solution will not be significantly modified 
by the self-attraction and only the fluctuations around it will be affected. This means that we can write down the 
(linearized) loop energy around the solution $ = 0, (j> = 4>i 00 p in a way similar to the last section, namely 

/L/2A I-L/2X 
60T,|(5(M + ^F / 6<&T15<&dt + V (D c (<5tf)) (92) 
-L/2X J-L/2X 

The last term V (D c ) that we introduced here - in accordance with above stated assumptions - represents the 
interaction potential of two overcrossing parts of DNA that have a closest distance D c . To keep the problem tractable 
we approximate here the distance D c by the perpendicular distance of the two crossing DNA parts at the equilibrium 
( mean ) crossing point t c of the homoclinic loop 

/t c />t c 
sin 5$ (t)dt w A / 8<&{t)dt (93) 
-t c J— to 

The crossing point t c will be given by the (in-plane) projected self-crossing of the loop. This implies the condition that 
the integral (over the interval [— t c , t c ]) of the x-component of the loop tangent vanishes, i.e., J_° cos<j)i oop (t) dt = 
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which leads to the following implicit equation for t c : 

f— E(t c /y/m\m)) 

t c = Vm (94) 

1 — TO 

Before we compute further it is interesting to have a short look at D c from Eq. ^Sl Because D c depends only on 
the out-of-plane perturbations, 6$, the in-plane (5(f>) problem stays unaffected. Note further that the out-of-plane 
rotational mode do (the generator of an infinitesimal rotation) leaves the distance D c unaffected: formally because 
$o (t) is an odd function, physically because rotations leave distances fixed. 

Now the partition function resulting from Eg. 1921 for any given V can be written as follows 

Q v = Q 2D Ql (95) 

Only the out-of-plane partition function is modified by the presence of the contact potential and is given by 

poo 

Ql= / e-m*JQ ± (# c )d0 c (96) 



where we introduce the angle § c — D c /2t c X that measures the (over the loop) averaged angle deviating from the plane 
i? = 0. The expression Q± (i? c ) denotes the properly constrained partition function 

Qx (tfc) = [*(#c-^- f C 60dt) e " i ^r Z f-L% «**is**t v [M ] ( 97) 



J \ 2i c J _ tc J 
To compute this path integral we replace the 8 function by its Fourier representation 



S U c _ _L J c 5 0dt) =^l°° e tp ^~^ IX «*) dp (98) 



The integral in the exponent is more elegantly written as a scalar product of 8$ with a " boxcar" function n (t) = 
H (t + t c ) - H (t - t c ) with (H (x) = 1 for x > 0, H (x) = 1/2 for x = and H (x) = otherwise): 

t c j-L/2X 

8-ddt = IL(t) 8d (t) dt = (n|&?) (99) 

-t c J-L/2X 

where we introduced the scalar product (f\g) = /f^A / (*) 9 CO ^ ^ n n °t a tion and by virtue of Eas. I98landl99l 
the partition function Q± (D c ), Eq. can be recast in a more transparent form: 

Q± (# D ) = J_ /°° e iP*. /'e-^< w l*- L l w >- i *< n l w >2?[<ytf]dp (100) 
2tt ./-oo 7 

We have now to compute the following path integral 

Ql(p) = [ e-^w^w-^w^vM (101) 



This path integral can be rewritten in as a Gaussian path integral in the presence of an external source current 
j(p, t) = ipLJ(t) coupled linearly to the fluctuating field 8§. We refer to Appendix B for the computation of this kind 
of path integral. The important point is that it can be written in the form 



^■ 7r - u ± \ 2A ' 2A 



where the determinant Dj_ (— jt, ^r) is independent of the source term. Along similar lines as in Eq.|SUwe go to the 
limit k — > 0: 

Q ± (p) k = q (j~\™)fj 4*W<k) Q ± (p) (103) 

where the Jacobian is given by expression Eq. 1701 



20 



It is very difficult task to calculate Eg. 11021 for the case of a finite chain length, so that we restrict ourself to the 
limit of very long chains. In this case the determinant is given by Eg. IB 121 and the functional E[j] by Eo. lfj20l that is 

3 (3i* - 10) 2 . . 

PElj] = -— V ; V 104 

The implicit condition on t c , Eg. 1941 becomes in this limiting case 

t c = 2tanh< c (105) 

that has t c w 1.915 as the numeric solution. This corresponds to t he actual loop circumference of 2 x 1.915A. The 
Jacobian Eq. [70] in this limiting case obeys J _1 (m = 1) = 2-^/2/3. Using the fact that the zero eigenvalue can be 
written as [1q (k) = y = (c — 1) (c + 1) we can compute the partition function Eq. 11031 in the infinite long molecule 
limit for k — > (i.e. c = 1): 

g i( p) = l8 ^ e -^ e ^(^-io)p 2 (106) 
Note that because of the unstable mode Q± is imaginary. Transforming back into real space yields 

Qx (0c) = ^ J°° e^Q^ (p) dp = 8^ (l^j V e~& exp (f l f (tf c ) 2 ) (107) 
where we introduced the scale-independent (negative) elasticity constant for the out-of-plane tilting 

f = 373^10) ^ 108 ) 
Using Eg. 1961 we can deduce for any given (reasonable) potential V (z) the out-of-plane partition function 



Ql=«/r(iA 3/2 e -± r e rif w-m*) dz ( io9) 

A V 7r \ A / J_ 00 

where we introduced T — T /At 2 c and the out-of-plane distance z = D c = 2t c \d c . This result indicates that the larger 
the perpendicular distance D c the larger the partition function. This is intuitively clear as the system without the 
constraint " D c = const." is intrinsically unstable tending to increase the distance D c . We can now deduce the full 
partition function Q v by combining Eos. 1521 and fl 091 

Q v = ^0zjF fFL-**-i J 00 f'W-firWte (110) 

This expression has to be taken seriously only for sufficiently fast growing interaction potentials V (z) for which the 
integral above stays finite. Otherwise the system is metastable and the integral diverges. But even in the case when 
the bound state, say z — zq, is just a local metastable state the integral above Q v still makes some sense if the V (z) 
is very deep. In this case the system can be considered as being in quasi-equilibrium (on some experimentally relevant 
timescale). If for instance we approximate V (z) locally by a quadratic potential V (z) — (z — zo) we obtain 

^ v 32V2TLd /2 SMl£ x i BFT H ' P l 32V2TWp 2 E1p z ? BFr H 'r L > s 

The last expression is valid in the limit of strong localization, i.e., for A'A 2 V AF. 

Let us finally write down the force-extension relation resulting from the general expression Eq. 11101 We express the 
free energy G v = — /3 _1 \nQ v in terms of F and A (instead of A = ^/ A/F): 

flG- = -fFL + (L + 8!,, Jl- In (I t~ ^^l*) - 1„ ( ^i^f'^\ (112) 



This leads to the force extension relation 



L 2\ L J p^AF 4f3FL 2 L\ f^ e P rA - 1/2F3/2z2 -P v ^dz [ ' 
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The last term ~ (D 2 C ) / (LX) is negligibly small because (-D 2 ) scales typically as the squared polymer cross-section 
(for a short ranged surface contact interaction). In the most extreme case the contact distance D c could become 
comparable to A (the loop head size), i.e., < A 2 but the latter is still much smaller than LX. That means that 

(for reasonable parameters of F and L) the force extension relation of a DNA loop with attractive contact interaction 
will essentially be independent of the concrete realization of the self-interaction potential V (x) and we recover the 

result of the previous section Eq. |^ The last term O ^(/3FL) -1 ^ is always negligible for large forces and therefore 

the average extension is to a good approximation given by the first two terms of Eq. 11131 



The second term is the usual "straight WLC" fluctuation contribution in 3D, Eq. ^1 the last term is the force 
extension signature of the DNA loop, cf. also Eq. 1551 

This computation shows that force-extension relation is fairly independent of many details like how we stabilize 
the loop in 3D. The physical reason for this simple decomposition has its roots in the fact that the WLC fluctuations 
(leading to the second term in Ea. lll4|) and its state of deformation (third term in Ea. lll4f) couple only negligibly in the 

large force regime (giving merely rise to weak logarithmic corrections in the free energy and negligible O ((/3FL) ^ 

corrections in the extension (Az)). 

The force-extension relation, Eq. 11141 has the same functional form as the usual WLC expression, Eq.^J but with 
an apparent persistence length 

(l + 8if) 

which is Eq. ^ °f the introduction. This shows that one has to be cautious when one probes the stiffness of a stiff 
chain via a stretching experiment: if the chain contains a loop then one will infer from the data a value for the chain 
stiffness that is too small. This is obviously mainly a problem in cases when the contour length of the chain is not 
much larger than its persistence length. But even for L/lp = 10 one finds l a p pp w 0.3Hp and for L/lp = 50 there is 
still a remarkable effect, namely l p pp w 0.74^p. 



V. CONCLUSION 



We have calculated the partition function of DNA under tension featuring a sliding loop via a path integration in 
the semiclassical limit (i.e., on the level of a saddle point approximation). This path integral can be mapped onto the 
QM harmonic oscillator with a time-dependent frequency In this analogy the time-dependence reflects the shape of 
the DNA chain. As it turns out the planar ground state solutions (Euler elastica) are always just "simple" enough to 
allow the exact solution of the corresponding path integral. The special choice of the parametrization of the tangent 
vector to the DNA has made the application of the semi-classical approximation possible as the singular measure term 
(due to inextensibility constraint) has been found to be negligible in this case. 

Within the semiclassical approximation the equation of state of looped DNA under tension for very stiff polymers 
is valid for any value of the applied force. The force-extension relation has been found to be expressed in terms of 
Jacobi elliptic functions and the force-extension curve provides two different scalings for weak and strong forces. For 
long DNA chains, the semiclassical approximation is valid only in the regime of strong stretching. In this force regime 
we proved that the elastic response of DNA is (up to logarithmic corrections) indistinguishable from the response of 
a non-looped WLC with the same contour length but a smaller persistence length. As we demonstrated the entropic 
fluctuations of the system are only marginally affected by the DNA shape, i.e., the entropies of the overall straight 
and of the looped conformation are essentially the same. What changes considerably when going to the looped state 
is the enthalpic part. It is the latter contribution that causes the apparent renormalization of the chain stiffness. This 
remarkable effect suggests that the results of corresponding micromanipulation experiments have to be interpreted 
carefully, especially in the case when the contour length of the chain is on the order of its persistence length. 

The looped DNA chain that we presented here should be considered as a paradigmatic model case. We believe that 
in the future this powerful approach will be applicable to a wide range a problems regarding semiflexible polymers. In 
Ref. |23| we already applied this method to DNA chains bearing deflection defects. Analytical results were obtained 
in the large force limit for experimentally interesting situations, e.g. for DNA with a kink-inducing bound protein 
and the problem of anchoring deflections in the AFM stretching of semiflexible polymers. Expressions relating the 
force-extension curve to the underlying loop/boundary deflection geometry were provided and applied to the case of 
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the GalR- loop complex |24j. The theoretical predictions were complemented and quantitatively confirmed by MD 
simulations [23j • Another non-trivial application of the semiclassical formalism concerns the buckling of rigid chains 
|48| . e.g. of microtubuli 



APPENDIX A: THE OUT-OF-PLANE DETERMINANT 



In this appendix we compute the out-of-plane determinant in the presence of a small magnetic field 
Dj_(— L/2X, L/2X) associated to the following out-of-plane fluctuation operator (cf. Eg. I66|) 

d 2 2 / * ■ N \ f 4 + m\ k 



In order to obtain a usual Lame equation we consider the transformation t' — > so that Tj_ can be written as 

1 / d 2 



rn 



T l = — ( —^2 + 6msn2 W\ m ) - (4 + m) + 3fc ) (A2) 



dt' 2 



Here we introduce 3k /m = k/F for later convenience. It can be shown 41] that the determinant is given by a 
particular solution of the following generalized second order Lame differential equation 

Tly(t') = (A3) 

Specifically, the solution satisfying the following boundary conditions: 

dy 

y(-K(m)) = 0, and — \- K(m) = y/m (A4) 

gives the desired result via f x( — 1X> = y{K{m)). 

We now solve Eq. IA3I with a method suggested in Whittaker and Watson's book ^3] • We first rewrite Eq. IA3I in 
the form 

y" - (6msn 2 (t\ m) - e) y = (A5) 
with e = 4 + m — 3k, and introduce the periodic variable 

z = cn 2 (i|m) (A6) 

in terms of which Eci. I A5I becomes 

p (z) y" +q(z)y' + r{z)y = (A7) 

with p (z) = —mz 3 + (2m — 1) z 2 + (1 — m) z. Now consider two linear independent solutions y\ (z) and ?/2 [z) from 
which we build the function M (z) = y\ (z) yi (z). One can then prove that this function satisfies the following third 
order differential equation 

2p (z) M'" (z) + 3p (z) M" (z) + {p" (z) + 2 (e + 6m (z - 1))) M' (z) + 6mA/ (z) = (A8) 

whose solution is a simple periodic function of the form M (z) = z 2 + az + b with the coefficients 

4-2m-e , , (1 + m - e) (4 + m - e) 

a = and b = i '- (A9) 

dm 9m z 

Because the Wronskian W of Ea. lA7l is given by W = y\y<i — yiVi = , C we can deduce that the solutions 

Wz(l — z)(l—mz) 

of Eq. IA7I are necessarily of the form 

yi 2 (z) = JITTzjexp I — j | (A10) 



with 



C = — ^ [(4 + m - e) (1 + m - e) (1 + 4m - e) (e 2 - 4 (1 + m) e + 12m)] 1/2 (All) 
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Introducing the transformation z = 1 — sin <p the integral in Eq. IA10I can be rewritten 

7=2 [ - # 

J v 1 — m sin 2 cf>M (sin 2 (ft) 

With the help of the following fractional decomposition 

1 Di D 2 



M (sin 2 



1 



1 



where the coefficients are given by 



2 + a - \Ja 2 - Aab 



Bo = 



2 + a + Va 2 - Aab 



and 



Di = 



2 + a- Va 2 - Aab 



D 2 = 



2 + a + Va 2 - Aab 



2{l + a + b) ^a 2 -Aab 2{l + a + b) Va 2 - Aab 

Eg. IA12l can be written in terms of the Jacobi elliptic function of the third kind II [n; (p\m) 



J{z) = 2L>iII 



Si 



; arcsin y/T— z\r 



2£) 2 n 



Bo 



; arcsin yl — z\m 



with 



II [n; ip\m] 



(l — n sin 2 (ft) (l — m 



sin 



,V2 



Then formally the solution of Eq. IA5I can be written 



J/1,2 (t) = y/M (t) 



exp 



(— 



J {t) 



(A12) 



(A13) 



(A14) 



(A15) 



(A16) 



(A17) 



(A18) 



At this point we mention that for n > 0, C is complex. Then the solution satisfying the boundary conditions Eq . IA4I 
is given by the following linear combination of two solutions given in Eq. IA10I 



J(-K)-J(t) 



(A19) 



This solution is valid only the interval —K [m] < t < because of relation Eci. I A6I between z and t. Note that 



J(-K[m]) = 2DJA 











+ 2L> 2 n 









(A20) 



In order to compute the determinant we need the solution y (t) for < t < K [m] that is also a linear combination of 
solution Eq. IA10I 



y(t) = -Vm 



\C\ 



^ainffl(J(i) + J(Jf[ m ])J 



Then for t = K [m] we deduce the determinant Dj_ = y (K [m]) so that 

-M{K[m\) . 



\C\ 



sm(\C\J(K [m])) 



(A21) 



(A22) 



This expression is valid for any value of k and could be used for the study of looped DNA in strong magnetic fields. 
Here instead we are interested in considering the limit of very small k (or, equivalently, k). For this we consider the 
expansion of the Jacobi elliptic function. 

Expansion of II ^-\m : 



24 



[li: 



From Ea. HTl we deduce the expansion B x w -^^k + O (k 2 ). B x being negative we find the following relation 

n (m — 1) 



n 



n = — — , m 
Si 



II[A,to] + K[m] 

(1 — nj (m — n) m — n 



(A23) 



with TV = yr^ 1 . As m < N < 1 we have also IT [N, to] = K [m] + ~5 2 (1 - A (e)) where A (e) is Heuman's Lambda 
function with 



£ = arcsm , 



(1 — to) 



1-N ■ / 1 1 n rr 

= arcsm \ j - w — y (1 — m) k 



1 — n to 



and 



N 



(to — n) (n — 1) 



(1 — N) (N — to) y n(l-TO) 2 (1 - to) ^(l - to) Jfe 



(A24) 



(A25) 



^From the relation Ao (e) = ^ {AT [to] E [e, 1 — to] — (A [to] — A [to]) A [e, 1 — to]} and using the two expansions 
E fe, 1 — toI w £ and A [e, 1 — toI ~ e, we deduce 



Ao (e) « -£ A [to] » ^^(1-™)^ 

7T 7TTO 



(A26) 



Now Eqs. IHT1 abd EH we the expansions 



which allows to write finally 



D 1 



2|C|Ain 



and \C\a-y/k(l-m) 

(1 — TO) fc TO 



A, 



2 J (1 TO)fc (A [to] - A [to]) - 7T 



(A27) 



(A28) 



Expansion of II j^\m 

^From Eg. IA14l we deduce the expansion B 2 ~ 1 + Since B 2 > 1 and to < -^A < 1 we can use the relation 

1 



n 



n = — , to 
Bo' 



= A [to] + -8 2 (1 - Ao (e)) 



where 



1-n 1 



e = arcsm , 



(1 — to) to y (1 — to) 



and 



so that finally 



62 = 



(l-ra)(n-m) ^(l -m)k 



\ / \ 2 , 2A[to] / fc 
A £ « -£ A [to] « — -W 7: r 

TT 7TTO 1/ (1 — TO) 



As A 2 « -to 2 + O(fc), we find 



2|C|A 2 n 



' 1 V 




— \m 


,1./ 




to y 



(A [to] - (1 — to) A [to]) — TT 



(A29) 



(A30) 



(A31) 



(A32) 



(A33) 



25 



Collecting the two expressions Eg. IA28| and Eg. IA33l we deduce that 



sin(|C| J(K [m])) w — W ., ((2-m)E[m] - 2 (1 - to) K [to] - 2tt) (A34) 

to y (1 — to) 

With this result in hand and with the expansion M (2K [to]) w -K we deduce that the determinant Eg . IA22I becomes 
in the small k limit: 



which writes in terms of k 



D1 » —^rr (2 (l-m)K [to] - (2 - m) J5 [to]) (A35) 

777/ (1 — 777) 

^ * 1 ^(1- m) (2 (1 " m) * H " (2 " m) ^ H) (A36) 

APPENDIX B: COMPUTATION OF THE PERPENDICULARLY CONSTRAINED PARTITION 
FUNCTION IN THE LIMIT OF INFINITE LONG CHAIN 

Here we evaluate the path- integral given by Eq. I1U1I It is equivalent to a special realization of the path integral of 
a QM harmonic oscillator with a time dependent frequency u> (r) and a driving force j (t) : 

/ K j] = f (Xl ' Tl) V [x] = ( (Xl ' Tl) V [x] e *Ctt W^-WW)^/,? ^M^r) (B1) 

The latter can be computed exactly (cf. Refs. p^l 



V 2mHD (n,To) 

The first factor on the rhs of Eo. IB2I represents the fluctuation contribution. Here D{t±,to) is the functional de- 
terminant of the (j-independent) operator T = cP/dr 2 + u 2 (t) normalized by the free-particle operator cP/dr 2 



n - t \ d 2 /d 

The second term in Eo. IB 21 involves the classical action S[j, x c i] where the j-dependent classical path x c \ (t) is the 
solution of the corresponding Euler-Lagrange equation 

motel (t) + mu 2 (t) x cl (r) = j (r) (B4) 

with boundary conditions x c i {tq/i) = a^o/i- Using Eo. IB4l the classical action can be rewritten as 

777 1 

s a, x c i] = y xci ( t ) ±ci ( t ) 1^0+2 y j ' ^ Sd dr ( bs ) 

Now in our concrete case we have to evaluate 

Q ± (p) = / P [M] wruwt+Ji 1 i(t)«* (B6 ) 

7(0, t ) 

with (asymptotic) boundary conditions at i]/ — ±00. The operator T* is again given by 

T ^VdT 2 -^^W + c ) (B7) 



20 



with c = yjl + y and N = 2. The source term is given by j (t) = — (H (t + t c ) — H (t — t c )). In our case the inte- 
gral IB2I fafter "Wick rotation" t — > — it and the replacement l/h — ► (3, m — > V 'AF, w 2 (t) — > (l — 6/cosh 2 (t) + k/F) 
etc.) has the following form 



with the classical action given by the following expression 

PS \j, W d ] = (t) d (t) |^ - i r j (r) tf d (r) dr = -f= t t G (t, if) dtdt' (B9) 

1 1 J-oc 2p\J At J-t c J-t c 

We first compute the fluctuation determinant D (— nr,jx) via the Gelfand-Yaglom method. It states that 
D 2x) = / (*) lt=i/2A where / (<) is the solution to T^/ (t) = with initial values / (-^) = and / (-^) = 1. 

The solution can be written in terms of the two linearly independent solutions (cf. |45j^l 

/I A \ W+1 
,JV+1 u\\ 1 „±ci 



/ ± (,) = cosh- M ^___J e - (BIO) 

For N = 2 the two independent solutions write 

/i (t) = e ct ((c - 2tanht)(c - tanht) - cosh -2 (*)) 

/ 2 (t) = e~ ct ( (c + 2 tanht)(c + tanht)- cosh -2 (*)) (Bll) 
and the general solution is given by / (t) — C±f (t) 1 + C%f (t) 2 . The Gelfand-Yaglom initial conditions in the limit 

L_ 

2\ 



L 1 (where we may safely set tanh (±<rr) ~ ±1, cosh 2 ss 0) determine the coefficients C\ and C2 



Ci - — ; %i ^ 

2c(c + 2)(c+l) 

2 2c(c-2)(c-l) 
Evaluating / (t) at the right boundary t — -kr we obtain for L/2A 3> 1 

x V 2A'2Ay 2c(c + 2)(c+l) V ' 

Note that unlike for the in-plane operator Qn case (where a close to eigenmode appears and creates artifacts) here 
we need not to renormalize D 1 ^ (— ttt, jrj an d Q± as long as the value of c is larger than 2. 

To compute the classical action Eq. IB9I consider the Euler-Lagrange equation which reads here 

pVAFT K J§ cl {t)=j{t) (B13) 

with the boundary conditions 5"d c i (±00) = 0. To solve this inhomogeneous differential equation we construct the 
Green's function [46| G{t,t') that is the solution to 

TIG(M') =S(t-t') (B14) 

with G (t, t') — G (t' , t) and proper boundary conditions G (±00, t') = 0. The latter gives the solution to Eg. IB 131 via 
the simple convolution 

Wd (*) = (pVAf) 1 f G (i, t') j (*') dt (B15) 

J — GO 

For our Dirichlet boundary conditions the Green's function generally writes 

q (fj tl) = H (f - t) h (0 h (t) + H(t- f) fi {?) h (t) (B16) 
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with /i two /2 being two (arbitrary) linearly independent solutions to the homogeneous equation Tj_f = satisfying 
the (one sided) boundary conditions /i (— oo) = and /2 (oo) = respectively. The constant W is the Wronski 
determinant of the two solutions, i.e. 

W = h (*) h (*) - fx (*) h (*) = const. (B17) 

We already know the two solutions, cf. Eq. IB11I Their Wronskian IB 171 is given after short computation by W = 
—2c (c 2 — 2) (c 2 — l). Inserting that and Eq. IB 111 into Eq. IB 161 gives a lengthy expression for G (t, t'). Fortunately 
there is no need for writing out explicitly neither G (t, t') nor <W C ; (t) as we are only interested in S [j, x c i] from Ea. lB5l 
(with x c i = Sd c i). Using Ea. IB 151 together with the boundary condition 5d c i (±oo) = leads to 

p 2 rU rtc 



pS[jM= A J f — / / G(t,t')dtdt' (B18) 
kpt%\/ At J-t c J-t c 

Inserting the Green's function, Eq. IB16I and exploiting fa (t) = fi (— t) we find 

PS [j > ^ = 4*^- 2) (<*-!) U H (0 \l h H) "'J M = 4^V-2) (ca -l) f (C ' tc) (B19) 



The involved double integral / (c, t c ) depends on the variable c = yT + y and the numerical constant t c in a compli- 
cated manner. But here we only need the case k — > 0, i.e., c — > 1. The expansion of the integrand around c = 1 (to 
lowest order) followed by the double integration gives (up to terms on the order of (c — 1)") 

I(c,t c ) = -6(c- 1) [2t c + 3 — % 5tanh£ c ) = (3t 2 c - 10) (c- 1) w 2.88 (c - 1) 

V cosh t c J 2 

Here we made use of the definition of t c , Eq. I1U5I The limit c — * 1 (i.e. k — > 0) can now be performed safely in 
Eq. IB19I and the action S [j, Stfrf} writes 

hm PS [J, Stfd] = _lM_B p2 ra _^2L p2 (B 20) 

This leads finally to 



w - ^ Ww^--*' h ^° y (B21) 
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